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A SIMPLIFIED MATHEMATICAL MODEL FOR CALCULATING 
AERODYNAMIC LOADING AND DOWNWASH FOR MIDWING 
WING-FUSELAGE COMBINATIONS WITH WINGS OF 
ARBITRARY PLAN FORM 

By Martin Zlotnick and Samuel W. Robinson, Jr. 

SUMMARY 


It is shown that, for the purpose of calculating aerodynamic loads 
on the fuselage, the midwing wing-fuselage combination in subsonic flow 
can be represented by a simple system of horseshoe vortices located on 
the wing with images located inside the fuselage. Using this simplified 
mathematical model, a method for calculating the lift and longitudinal 
center of pressure on the fuselage in the presence of the wing is 
presented. 

In addition it is shown how the simplified mathematical model can 
be used for calculating the downwash behind the wing and for calculating 
the spanwise lift distribution on the wing in the presence of the 
fuselage. 


INTRODUCTION 


Mutual interference between wing and fuselage has a significant 
effect on the pitching moment of the wing- fuselage combination since 
the longitudinal distribution of the aerodynamic loading on the fuse- 
lage is altered by the presence of the wing. For conf igurations with 
unswept wings Multhopp (ref. l) has developed a theoretical method for 
calculating the pitching moment on wing-fuselage combinations which gives 
good agreement with experimental results. However, for calculating the 
pitching moment on swept-wing configurations, only semiempirical methods 
such as that of reference 2 are available. 

The aerodynamic loading on the fuselage in subsonic flow is consid- 
ered in references 1 and 2 to be made up of two parts. The first part, 
which is due to the fuselage angle of attack resulting from both the fuse- 
lage geometric angle of attack and the upwash angle induced by the wing. 
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has been calculated in reference 1 for straight-wing configurations and 
modified in reference 2 for the swept-wing case. The second part, which 
is often referred to as the wing lift "carried over" by the fuselage 
(which will be referred to in this paper as the induced lift), has been 
calculated analytically for the straight-wing case by Multhopp. However, 
for swept wings Muithopp 1 s calculation cannot be applied and no other 
theoretical method has been available. In reference 2 this component 
of loading is estimated by an empirical method. 

In order to calculate the induced lift on the fuselage in combina- 
tion with a wing of arbitrary plan form, a method is suggested in the 
present paper which is based on Lennertz's (ref. 3 ) theoretical work. 
Lennertz's results, which are concerned only with an unswept lifting 
line passing through the axis of ah infinitely long cylindrical fuse- 
lage, are, in effect, generalized so that in subsonic flow the magnitude 
and center of pressure of the induced lift on the fuselage when combined 
with a swept lifting line, or even a lifting surface, can be calculated. 
A numerical example is given in the appendix to illustrate the method. 
The effect of finite fuselage length is estimated qualitatively from 
results of an approximate calculation of the variation with slenderness 
ratio of the induced lift on an ellipsoid of revolution combined with an 
infinite vortex. The approximate calculation is shown to give results 
which agree with results of Vandrey (ref. 4 ). 

In calculating the induced lift on the infinitely long fuselage, 
the wing-fuselage combination is replaced by a simple system of horse- 
shoe vortices (or doublets) on the wing, with images inside the fuselage. 
This representation may also be used, with some modifications, in calcu- 
lating the lift on the wing in the presence of the fuselage and also the 
corresponding downwash. However, for calculating the downwash and the 
lift on the wing, the simplified representation of the wing-fuselage 
combination, although it is considered to be adequate, is no longer 
rigorous, and the general applicability will depend on experimental 
verification. 


SYMBOIS 


A 

a 

b 

c 

c 


aspect ratio (wing alone) 

maximum radius of body of revolution 

wing span 

wing chord 

mean chord, S/b 
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c^ local lift coefficient 

d major axis of ellipsoid of revolution 

Lf total lift on the fuselage 

p pressure 

q dynamic pressure 

S wing area (wing alone) 

s semispan of hound leg of horseshoe vortex 

Uoo free-stream velocity 

Umax maximum longitudinal velocity on surface of body 

x C p longitudinal center of pressure 

r b vortex strength' 

7 spanwise loading coefficient on wing, ccj/c 

p mass density 

Subscript: 

n particular pair of horseshoe vortices 

A prime denotes that the quantity is dimensionless with respect 
to a; an asterisk denotes that the quantity is dimensionless with 
respect to b/2; and a bar over a symbol denotes that the quantity is 
dimensionless with respect to s. 


BASIC CONSIDERATIONS 

Midwing Configurations With Axisymmetrical Fuselages 


In the major part of the analysis the fuselage is assumed to be an 
infinite cylinder. To obtain a qualitative estimate of the effect of 
the finite fuselage length, the variation with slenderness ratio of the 
lift on an ellipsoid of revolution with a wing of infinite span will be 
calculated approximately in one of the following sections. An approxi- 
mate method for making the small correction for this effect will be 
indicated. 
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Review of results obtained by Lennertz . - Lennertz (ref. 3 ) has cal 
culated the lateral and longitudinal lift distribution on the fuselage 
of an idealized wing-fuselage configuration, in which the fuselage was 
represented by an infinite circular cylinder and the wing by a vortex 
having constant spanwise circulation. The vortices trailing from the 
wing tips have images inside the cylinder, and the bound vortex is 
extended inside the cylinder to join the trailing image vortices as 
shown in figure 1. With this conf iguration, the boundary condition on 
the cylinder of zero velocity normal to the surface is satisfied only, 
at infinity and in the plane normal to the cylinder axis which passes 
through the vortex so that it is necessary to superimpose an additional 
potential, which Lennertz calculated. For this case the lateral lift 
distribution .was obtained by considering the momentum change in a verti 
cal plane infinitely far behind the wing, and the longitudinal lift d-is 
tribution was obtained by the use of Bernoulli’s equation. 

Induced lift on an infinite cylinder with a finite wing having 
constant spanwise circulation distribution .- It will be shown that 
the fuselage lift Lf and its lateral distribution dLf/dy are not 
affected by the additional potential so that only the components of Lf 
and dLf/dy due to the vortex potential need.be calculated. 

Neglecting all other singularities, the lift on a longitudinal 
section of the cylinder dy (see fig. l) due to the vortex potential 
is calculated as follows: 

The pressure p of any point on the surface of the cylinder can 
be written 


P = P 0 + [(^00 + Aii ) 2 + Av 2 + Aw^j 

where Au, Av, and Aw are, respectively, the longitudinal, lateral, 
and vertical components of velocity induced by the vortices. The sec- 
tion lift dLf/dy is then written 


dLf • pco 

— - (p u - p.)a* 


where p u and pj are the pressures on the upper and lower surfaces, 
respectively. 
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Then 


dLf 

dy 


= 2pU 


Au dx 


( 1 ) 


Only the velocities induced by the bound Vortex contribute to the lift. 

It may be noted that, since the distribution of the longitudinal 
Velocities induced by the bound vortex is symmetrical about its axis at 
every section dy, the longitudinal center of pressure of the lift due 
to the vortex potential is on the axis of the bound vortex. ' 

% 

A closed expression for the integral in equation (l) may be readily 
derived as follows: Consider the rectangular path indicated by the 

dashed line in the upper right sketch in figure 1 and the cross section 
downstream at infinity shown in the lower right sketch of the same fig- 
ure. The line integral of the tangential velocity component, taken 

around the complete rectangle, ^Au ds must be where the path 

links one of the horseshoe vortices (|w| < |y l < l a 0’ and must be zero 
where' it does not link one of the horseshoe vortices (0 < lyl < p£_ ]. 

V 1 1 |b/2 / 

Then the desired integral in equation (l), which is the longitudinal 

portion of the complete line integral Au ds, must equal ^ Au ds 

minus the line integral along the short vertical line at infinity 

/ Aw(oo,y,z)dz, where u refers to the upper surface and l refers to 
'•'u 

the lower surface. The integral J ' Au dx is 


p 00 p pi 

/ Au dx = (^Au ds - / Av(oo ,y,z)dz 
J -oo d J u 


( 2 ) 


Since 


30 

Aw = - —— 

3 z 
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/ Aw dz = 0 U - 0 2 

where 0 U and 0^ are the potentials on the upper and lower surfaces, 

2 

respectively, in the plane x = oo. When < lyl < la I, 

b/2 11 
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and 


(^Au da = 0 


Then, from equation (2) , the lateral lift distribution of the induced 
lift can be written 


(-a < y < a) (3) 


and after integration over y 

Lf - 2pUJb(a - (4) 

which are the same as the expressions obtained by Lennertz for dLp/dy 
and the total fuselage lift Lf. The components of Lf and dLf/dy 
due to the additional potential must therefore be zero. 

The longitudinal lift distribution on the cylinder calculated by 
Lennertz, which included the effect of the additional potential as well 
as the vortex potential, is slightly different from that which would be 
calculated due to the effect of the vortex / potential alone. However, 
both distributions are symmetrical about/ the axis of the bound vortex. 

The longitudinal center of pressure of the lift due to the vortex poten- 
tial must be at the bound-vortex axis because the longitudinal velocities 
induced on the surface of the cylinder by the bound vortex are the same 
in front and in back of the bound vortex as noted before. That the 
longitudinal center of pressure of the lift due to the vortex potential 
and the additional potential must be at the axis of the bound vortex as 
the Lennertz calculations show (although he does not say this explicitly) 
can be shown as follows: 

Superposition of infinite vortices canceling the semi-infinite 
trailing vortices of figure 2(a) as shown in figure 2(b) will not change 
the longitudinal lift distribution on the cylinder since no longitudinal 
velocities are induced. It is apparent that the system is the same as / 
before, with the direction of the trailing legs of the horseshoe vortices 


dT i j 

dy 


3? - e u ~ r t 


l - 


— tan' 
jt 


2, 

■1 2 


- y^ 


(!) - 
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reversed; therefore, the longitudinal distribution on the infinite 
cylinder must be symmetrical about the axis of the bound vortex, since, 
if it were not, the longitudinal lift distributions on the cylinders in 
figures 2(a) and 2(b) would be different. 

From the analysis in this section of the wing and cylinder combina- 
tion it can be seen that, for calculating the lift on the cylinder, its 
lateral distribution, and the longitudinal center of pressure, it is 
necessary to consider only the effect of the vortex potential. This 
lift due to the vortex potential will be referred to hereinafter as the 
"induced lift." 


Simplified representation of the wing-fuselage combination . - In 
this paper, the wing-fuselage combination will be represented by a system 
of discrete horseshoe vortices and images, so that configurations with 
wings of arbitrary plan form may be treated. It is necessary to super- 
impose two pairs of horseshoe vortices of the type shown in figures 1 
and 3 to obtain the vortex- image system of figure 4. The vortex- image 
system of figure 4 can be used to represent wings of arbitrary plan form 
in the manner shown in figure 5 by locating the bound vortices on the 
wing quarter-chord line or by distributing them over the wing surface in 
the same mann er. Since the discrete horseshoe vortex becomes a doublet 
line as the length of the bound leg approaches zero, it is also possible 
to represent the wing by a continuous distribution of doublets. 

Figure 3 shows the larger pair of vortices with span 2(y n + s) and 
strength +r near the smaller pair of span 2 (y n - s) and strength -r. 
When the bound legs of the smaller are moved to coincide with the bound 
legs of the larger, the net vortex strength along the section where they 
coincide is zero, and the remaining sections of the larger bound legs 
form the bound legs of the desired system (see fig. b ) . From equa- 
tions (3) and (4), the lateral distribution of the induced lift and the 
total induced lift for the vortex- image system shown in figure b are 
given by the following expressions: 



(¥•) * 
x ' n 


± 2 (yn - B )^ a 2 - y 2 + _ ± 2 (y n + s)^a 2 - y 2 


|tan” x — - — — tan 

(y n - s ) 2 - a 2 (y n + a) 2 - 


2 - a 2 


(5) 
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ALf 

qs" 



(6) 


is the loading coefficient on the wing at the station y n . 


is known that the longitudinal center- of pressure of the 
induced lift (and total lift) for each pair of horseshoe vortices and 
images is still located at the axis of the hound vortex, it is possible 
by superposition to calculate the pitching moment on the fuselage in 
the presence of a wing of arbitrary plan form if the lift distribution 
on the wing is known, although the complete longitudinal lift distribu- 
tion cannot be calculated unless the component of loading due to the 
additional potential is calculated. The total induced lift and its 
lateral lift distribution can also be calculated by superposition. The 
method is discussed in one of the following sections, and an illustra- 
tive example is given in the appendix. 


where 


(?! 


Since it 


Effect of Finite Fuselage Length on the Induced Lift 

To obtain an estimate of the error involved in the assumption of 
the foregoing analysis that the fuselage is infinitely long, an approxi- 
mate calculation will be made of the effect of slenderness ratio on the 
induced lift of an ellipsoid of revolution combined with an infinite 
vortex. The limiting case of the spherical fuselage is treated first 
in the manner of Vandrey (ref. 4 ), and then the general case is treated. 

Calculation of the lift on a sphere-vortex combination .- The poten- 
tial of the sphere combined with an infinite vortex (fig. 6(a)) is written 


^ = — tan"-*- — + U x 
2« x 


1 + i 


(jx 2 + y 2 + z 2 ) 


( 7 ) 


where the potential for the sphere in the free stream is 
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®1 = 


('jx 2 + y 2 + z 2 ) 3 
and the potential for the vortex is 


w 


( 8 ) 


, , _i z 

Oo = — tan ^ — 
2 2n x 


( 9 ) 


The lift on a section of the sphere in the plane y = Constant is 


dLf f ^ 2 ' y2 / * 

5T ’ ( p " - p ') dx 


( 10 ) 


where p u and p i are the pressures on the upper and lower surfaces, 
respectively. From Bernoulli's equation. 


(fu + Au) 2 + (v + Av) 2 + (w + Aw) 2 


(11) 


and 


Pu - P i 


2p(u Au + w Aw) 


where Au, Av, and Aw are the velocities induced by the vortex on 
the upper half of the sphere and u, v, and w are the local velocities 
on the surface of the sphere. Thus 


Au = ~ z 


2 it n 2 


y 2 


Aw = — x 


2 it a 2 


- y 2 


u = | u 

2 to p 


w = Ju 22. 
2 00 2 
a 


( 12 ) 


RESTRICTED 



NACA RM L52J27a 


RESTRICTED 


11 


The lift on the section at y = Constant is 


dL f = 3 £b P U o» 
dy 2 * a 2 . y2 



z dx 


(13) 


which may be written in the form 


flT,p r f (y) 

— - = pUmax / 2 Au dx (l4) 

dy '-'-f(y) 


since 


Au 


£b z 
2n a 2 - y 2 


and since 


u = 


U, 


max 


when x = 0, so that 


U, 


max 


- lu 


For an infinite cylinder with an infinite vortex, equation (lit-) will 
also hold. In this case 


W = U M 


A r b 4 a 2 - y 2 
Au = — 


2it q 2 


y 2 + x 2 


and 


f(y) = 
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From equation ( 13 ) or ( 14) the lift obtained for the sphere is 


dLf 

dy 




and for the cylinder 


dLf 

dy 


= PUR 


co b 


Equation (l4) gives the exact result for the case of the sphere 
with an infinite vortex and the same result as the method of Lennertz 
for the case of an infinite cylinder with an infinite vortex. It is 
suggested, therefore, that equation (l4) may be used to obtain a quali- 
tative estimate for the intermediate case of an ellipsoid of revolution 
having a slenderness ratio between 1 and infinity. This assumption will 
be made in the following section and it will be shown that the value 
obtained for the induced lift by using equation (l 4 ) is very close to 
that obtained by the more accurate mathematical treatment of Vandrey 
(ref. 4) for the case of an ellipsoid having a slenderness ratio of 5 . 

Calculation of the lift on a combination of an ellipsoid with an 
infinite vortex .- From equation (14) the local lift on the ellipsoid ia 
written 


dLf 

dy 


= pU, 


max 


of(y) 
I- f(y) 


2 Au dx 


(15) 


From figure 6(b) , 

f(y) = ajl - (|) 2 


a f 

-© 

2 p 

-(■0 

rf 

2n 

x 2 + a J 

! i- (if - 

(if 


(16) 
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The integration indicated in equation (l4) yields 


1 ILf _ Umax/^oo 
p u oo r h d y 1 + | 


The values of Umax/^co given as a function of the slenderness ratio in 
reference 5 are shown in figure J. 

1 dLp z 

Figure 8 shows a plot of ■ — — against d/a. The result of 

a calculation by Vandrey for the case of an ellipsoid with §■ = 5 is 

shown to fall near the curve in figure 8. Since Vandrey’ s calculations 
appear to be accurate to 0.05, the agreement may be even better than is 
indicated in figure 8. Vandrey' s results are obtained by difficult 
computations for each body separately, and his method does not yield a 
general expression similar to equation (17). 


Figure 8 indicates that the induced lift on the ellipsoids having 
high slenderness ratios (d/a > 5) is about 90 percent of the induced 
lift on the infinite cylinder. Since the fuselage is similar effec- 
tively to a semi-infinite cylinder because of the wake which extends 
behind it, the loss in lift due to the finite fuselage length is about 
half of that indicated by the calculation for the ellipsoid of revolu- 
tion. The value of induced lift obtained by assuming the fuselage to 


be an infinite cylinder may be multiplied by the factor 



to correct for the finite fuselage length; however, this approximate 
correction can usually be neglected since it is nearly unity for most 
practical slenderness ratios. Since the correction is small, it is 
assumed that it may be applied directly for the finite-wing case without 


introducing significant error, so that the factor 

be multiplied by equations (5) and (6) to correct 
length. 




is to 


for finite fuselage 
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APPLICATION OF SIMPLIFIED MATHEMATICAL MODEL TO CAICUIATION 
OF THE AERODYNAMIC LOADING AND DOWNWASH 


In the following sections, the method for calculating the induced 
lift will he described and discussed and methods for calculating the 
spanwise loading on the wing and the downwash will be outlined. The 
methods for calculating the downwash and spanwise lift distribution 
are not rigorous since the effect of the additional potential, which 
has not yet been calculated, must be approximated by a simple correc- 
tion. Although the validity of the correction must depend on experi- 
mental verification, it is believed to be adequate, and the exact value 
of the additional potential may be incorporated into the method immedi- 
ately when it is calculated. 


Method for Calculating the Induced Lift on the Fuselage 

If the lift distribution on the wing in the presence of the fuse- 
lage is known, say from reference 6 or the method outlined in the 
following section (see fig. 9), the induced lift may be calculated very 
simply. A sample numerical calculation is shown in the appendix. It is 

only necessary to substitute into equations (5) and (6) the value of 

CC 7 

loading coefficient, — -, for each of the discrete horseshoe vortices. 

c 

The values of the increments of the lateral lift distributions A 

<iy 

on the fuselage due to all of the horseshoe vortices and their images 
obtained from equation (5) are superimposed to get the complete lateral 
lift distribution dLf/dy. The increment of total lift for each pair 
of horseshoe vortices and images AL^ acts at the bound-vortex axis, 
and the longitudinal center of pressure of the total lift Lf is 
obtained s imply by dividing the sum of the moments of the incremental 
total lifts ALf by the total lift Lf (the sum of the increments ALf ) . 

The results of the calculation for the lateral distribution and 
longitudinal center of pressure carried out in the appendix are pre- 
sented in figures 10 and 11, respectively. Figure 11 shows the incre- 
ment of total lift ALp contributed by each pair of vortices and images 
and it can be seen that the contribution to the total lift Lf due to 
the outboard part of the wing (the contributions farthest aft in the 
side view) are small compared to the contribution of the inboard part 
of the wing. The induced lift on the fuselage has also been calculated 
by using the wing-alone spanwise lift distribution shown in figure 9; 
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that is, the effect of the fuselage on the wing is neglected. It can 
he seen in figures 10 and 11 that, although there is about a 10-percent 
increase in the magnitude of dLf/dy and Lp due to the effect of the 

fuselage on the wing, the lateral and longitudinal load distribution on 
the fuselage is practically unaffected. The correction for finite length 
has not been included in the calculations. Its only effect would have 
been to decrease the magnitude of the lift, but it would not alter the 
lateral or longitudinal distribution. 

The total lift and moment on the fuselage can be obtained by adding 
the components of lift and moment due to the induced lift to the com- 
ponents of lift and moment due to the local angle of attack of the fuse- 
lage (which can be calculated as shown in ref. 2) . The lift and moment 
on the part of the wing outboard of the fuselage wall calculated by the 
method of reference 6 or the method of this paper (described in a later 
section) may then be added to the lift and moment on the fuselage to 
get the total lift and moment on the combination. 


Outline of Method for Calculating the Lift on the Wing in 
the Presence of the Fuselage 

The lift distribution on the wing will be calculated by equating 
the downwash angle induced at the three-quarter-chord line by the horse- 
shoe vortices centered on the quarter-chord line (see fig. 5) to the 
local angle of attack on the three-quarter-chord line at several points 
along the span. Since the boundary conditions on the fuselage are not 
completely satisfied by the vortex-image system of figure 5 in the 
region near the bound vortex, it is necessary to resort to certain 
approximations in calculating the downwash. Calculation of the exact 
values of these downwash functions would require the calculation of the 
additional potential which involves a great deal of time and effort, 
but if such calculations were made, the, values could be used directly 
in the present method and the restrictions suggested below would be 
eliminated. 

The approximations described in the following sections improve as 
the longitudinal distance from the downwash points to the bound vortex 
increases so that this method is considered to be best-suited for con- 
figurations having fairly small ratios of diameter to root chord, say 

for straight wings with _£iameter_-g 1 and for swep t W i n g S with 

Root chord d - 

Diameter < 1^ Calculations of the aerodynamic loading on a wing and 
Root chord 3 

tip-tank combination, using a method corresponding to the one described 
herein with the downwash points located about 1 tip-tank diameter behind 


RESTRICTED 



1 6 


RESTRICTED 


NACA RM L52J27a 


the bound vortex f % i), have been found to yield results in 

\Tip chord d l 

good general agreement with experimental results. 

Downwash in region near the wing bound vortex .- Although the 
boundary conditions are satisfied completely only at infinity, the 
trailing legs of a single pair of horseshoe vortices and images also 
satisfy the boundary conditions on the cylinder in the plane perpen- 
dicular to the cylinder axis which passes through the bound legs of the 
horseshoe vortices. In addition, the boundary conditions on the cylin- 
der are satisfied completely by the semi- infinite vortex-image system 
everywhere in the plane of the horseshoe vortex. Since the boundary 
conditions on the cylinder are satisfied exactly at the points noted, 
and are partly satisfied everywhere else, it will be assumed that the 
downwash due to the trailing vortices may be calculated approximately, 
at least in the plane of the horseshoe vortex, without introducing any 
correction factor. 

It will be necessary, however, to use a correction factor in calcu- 
lating the downwash due to the real and image bound vortices which have 
the greatest tendency to violate the boundary condition of zero velocity 
normal to the surface of the cylinder. It will be assumed that, along 
a line parallel to the bound vortex, the effect of the cylinder on the 
vertical flow induced by the bound vortex will be the same as its effect 
on a two-dimensional uniform rectilinear flow, so that the downwash 

induced by the bound vortices on that line is increased by the fac- 
a 2 

tor 1 + 

Ya 

The downwash angle at a point at y = y a and x = Xg (shown in 
fig. 4 ) is then written 


cl =j-iy 

^ 4«A 2s ^ 


r n 


n 


F nU,y,°) + G n(x,y,0) (l + sL\ 

V y a 1 


(18) 


where the downwash factor due to the trailing vortices F n and the 
downwash factor due to the bound vortices G n may be calculated by 
the Biot-Savart law. 
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The function F n is the sum of four terms, each having the form 


1 

*a - *n 

y n - y a + 1 

~ ya + x ) 2 + (*a “ %) 2 


1 


7n _ ^a 



+ 1 


(the term corresponding to the trailing legs of the horseshoe vortex 
centered at y = y n ) and the function G n is the sum of four terms, 
each having the form 


1 

y n - 

y a + 1 

y n 

- y a - 1 


x a - % 

^ - y a + 


1 

if>> 

1 

i) 2 . fc 

, - in) 2 


(the term corresponding to the hound leg of the horseshoe vortex 
centered at y = y n ) . 

Scheme for calculating spanwlse lift distribution on the wing . - The 
lift distribution on the wing is calculated by equating the downwash 
angle on the wing induced by the vortex- image system to the local angle 
of attack on the three-quarter -chord line of the wing at n points as 
in reference 6. Thus equation (l8) is written for each of the n sta- 
tions, and the simultaneous equations can then be solved for 7 n at 
n points on the span.: The + marks in figure 5 indicate the points 
where the downwash is equated to the local angle of attack. 

The effective angle of attack Oq is equal to the geometric angle 
of attack on the wing cig plus the angle of attack induced by the 
fuselage 

(% = a g .+ otf 1 (19) . 
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and 


of' = 



- Of 


( 20 ) 


where Of is the geometric angle of attack of the fuselage and the 


factor 1 + 


takes into account the increase in the vertical velocity 

of the free stream in the neighborhood of the fuselage as calculated by 
assuming the fuselage to be an infinite cylinder in a two-dimensional 
uniform rectilinear flow of magnitude Vctf. 


() The effect of finite fuselage length (referred to in ref. 6 as the 
inflow effect") must be included separately by multiplying 7 at each 
spanwise station by the factor 1+25 where 5 is the ratio of the 
local increment of longitudinal velocity due to the fuselage to the free- 
stream velocity. As shown in reference 6 , the factor 1+25 is used 
to account approximately" for the small increase in dynamic pressure in 
flow over the wing due to the increase in the local longitudinal 
velocity near the surface of a fuselage of finite length. 


Outline of Method for Calculating the Downwash 


The downwash calculations made in this section are for the case 
when the wing is at an angle of incidence with the fuselage at zero 
geometric angle of attack. The accuracy of the calculated results must 
be verified by experimental results; however, by comparison with results 
of downwash calculations for the wing alone, the results of the calcu- 
lations for the downwash in back of the wing-fuselage combination may 
serve to. give useful information regarding the effect of the wing-fuselage 
interference. 


rnmh-^ S +r S n0ted pr ® v f ousl y^ the "boundary conditions on the wing-fuselage 

■ the e wir the system 

,, , . ^ a 8 re ^t distance behind the wing, therefore the 

heoretical value of the downwash can be calculated accurately for a’ 
given spanwise lift distribution. .This calculation is simple^ and ?he 
results may be useful for certain applications. However the annroxi 
nations made m the preceding section for calculating the downwash at 
the wing three-quarter-chord line may give more accurate results for 

Se b“r?:x! Ulated ^ ^ Plane ° f thS Ving in the nearer 
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If the lift distribution on the wing is known, say from the method 
of reference 6, or the method outlined in this paper, the downwash angle 
in the plane x = » may be calculated by adding the downwash due to 
each pair of trailing vortices and their images. Thus, from the Biot- 
Savart Law and figure 4, the downwash angle at the point y = y , 

z = z a , and x = «. is given as 


Oa - 


2nA.s* 


z 

n 


?n F n 


(21) 


where 


F = 
r n 




(y n - + x ) + 


; 2 



y n + y a + 1 


(y n + y a + x ) 2 + 2 a 2 


y n + y a - 1 


(y n + y a 



+ 



- a^ 

y a + tr- + 1 

a y n 


y a 



1 



( 22 ) 


and 



is the local lift at station ±y n . 
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. The representation of the wake by discrete vortices is satisfactory 
if the downwash is calculated at points halfway between the trailing 
legs ( y a = y n ^. The values calculated at these points may be faired 

to obtain a continuous spanwise distribution of downwash angle. 


CONCLUDING REMARKS 


For the purpose of calculating the longitudinal loading on the fuse- 
lage in subsonic flow, a midwing wing-fuselage combination has been repre- 
sented by a system of discrete horseshoe vortices and images. Using this 
simplified mathematical model, a method for calculating the lift on the 
fuselage induced by the wing is presented and illustrated by a numerical 
example. This "induced lift"' can be added to the lift on the part of 
the wing outboard of the fuselage and the lift on the fuselage due to 
the upwash induced by the wing to get the total loading on the wing- 
fuselage combination with an axisymmetrical fuselage and a wing of 
arbitrary plan form. 

In addition to the method for calculating the induced lift, which 
is theoretically rigorous, methods for calculating the downwash far 
behind the wing and for calculating the spanwise lift distribution on 
the wing in the presence of the fuselage are outlined. In calculating 
the spanwise lift distribution on the wing, approximations are made to 
account for the effect of the "additional potential" so that the method 
is not rigorous. However, the effect of the additional potential may 
easily be incorporated into the method when it is calculated. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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APPENDIX 

ILLUSTRATIVE EXAMPLE 


A numerical example is given to illustrate the method presented in 
the body of the paper for calculating the induced lift on the fuselage. 


Geometric characteristics of the midwing configuration used in the 
illustrative example .- The plan view of the wing-fuselage, combination 
is essentially the same as shown in figure 5. The geometric data are 


Aspect ratio . 
Taper ratio 
Sweepback, deg 

a* 

s* 


8 

0.45 

45 

0.10 

0.05 


Spanwise loading on the wing in the presence of the fuselage .- The 
spanwise lift distribution on the wing in the presence of the fuselage 
tabulated below and shown in figure 9 is the same as the one obtained 
in reference 6 for a high midwing configuration. It. ia, however, quali- 
tatively correct for the pure midwing case and is used to illustrate 
the procedure for obtaining the induced lift on the fuselage. The lift 
distribution on the wing in the presence of the fuselage plotted in fig- 
ure 9 is tabulated as follows: 


n 

y n* 

x n* 

(?). 

i 

0.15 

0.15 

0.369 

2 

• 25 

.25 

.366 

3 

.35 

.35 

.356 

4 

.45 

.45 

.338 

5 

.55 

.55 

.319 

6 

.65 

.65 

.300 

7 

.75 

.75 

.266 

8 

.85 

.85 

.241 

9 

.95 

.95 

.200 


Longitudinal distribution of the induced lift .- Equation ( 6 ) which 
gives the lift on the fuselage due to a single pair of horseshoe vortices 
at and their images can be written 


i 
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The total lift Lf/qS is then 



The longitudinal center of pressure referred to the intersection of the 
quarter-chord line and the fuselage axis of symmetry (see fig. 5) is 



The computed values of Lf/qS and x^-p ' are obtained from the values 
presented in the following table: 
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= 2.51 


The values of 



are plotted as vectors In fig- 


ure 11. The location of x C p' is also shown. 


Lateral distribution of the induced lift . - Equation ( 5) which gives 
the lateral distribution of the induced lift on the fuselage due to a 
single pair of horseshoe vortices at (%'jtyn') and their images can be 
written 



The total lift at a station ±y' on the fuselage is then 



(A5) 
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Figure 9 *- Lateral lift distribution on wing alone and on a wing-fuselage 
combination. Sweepback angle, 4 5 °i taper ratio, 0.^5; A = 8; a* = 0.10 



Iculafed *✓///> fuselage effect on wing 
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